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Abstract. We present timing and spectral results of PSR B1937+21, the fastest known millisecond pulsar (P ~ 
1.56 ms), observed with RXTE. The pulse profile, detected up to ~20 keV, shows a double peak with the main 
component much stronger than the other. The peak phase separation is 0.526 ± 0.002 and the pulsed spectrum 
over the energy range 2-25 keV is well described by a power law with a photon index equal to 1.14±0.07. We find 
that the X-ray pulses are closely aligned in phase with the giant pulses observed in the radio band. This results 
suggest that giant radio pulses and X-ray pulses originate in the same region of the magnetosphere due to a high 
and fluctuating electron density that occasionally emits coherently in the radio band. The X-ray events, however, 
do not show any clustering in time indicating that no X-ray flares are produced. 
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1. Introduction modeled either with curved or broken power laws (Becker 

& Aschenbach 2002). 

Millisecond radio pulsars (MSPs) differ from ordinary pul- ^c,^ nmoTioi c *. Hon j- j m i 

r i • i . . , . , -, „„ PSR B1937+21 was the first MSP discovered (Backer 

k^l ' sars lor their short spin period m the ranee 1.5-30 ms , , -, n0 r,N , • , c -, „ . 

^ . . . . • "„_,„ ct al. 1982) and. with the period ol 1.56 ms, it remains the 

■ and tor their low spm-down rate of P < 10 J s s . , .,, ... , , , r™ 

* most rapidly rotating neutron star presently known, the 
CZ ' They are presumably very old objects, with spin down , ,. , , c , , ,. 

- - i J „q ,„in , , r • c i , distance estimated from the observed dispersion measure 

ages ol r ~ 10-10 ' yr and low surface magnetic holds m ,,> , r j i r ai. n i j.- r i A j- 

° „ , n _ . ... . , ,, , (DM) and from a model for the Galactic free electron dis- 

B ~ 10 -10 ' G. MSPs are believed to be recycled by hav- rm i o n j mno n j p t <m no\ 

, tribution (Taylor & Gordes 1993, Gordes & Lazio 2002) 

me accreted mass and angular momentum from an evolv- . „ „ , , , ,. AT . , , onno \ 

. . /.I i a, „ w r i, 1S 3.6 hpc (sec also discussion m Nicastro et al. 2003). 

me companion star (Alpar ct al. 1982). About 50% of all T , . , , . ., . a , , ln36 _i , 

, , v , „ Its spin down luminosity is K ~ 1.1 x 10 JD erg s and 

X-ray detected rotation-powered pulsars are MSPs and , ,. ~ , , , c 

J i , the dipolar magnetic held component at the star surface 

their X-ray emission seems to have characteristics similar . A1 1 x ^ , , n , , , 

J is ~ 4.1 x 10 G. Like the Grab pulsar (Lundgren et al. 

to ordinary pulsars with both thermal and non-thermal lgg5) pgR B0540 _ 6g ( Johnston & Romani 2003) and the 

components. The former can be described by a modified Qther Mgp pgR B1821 _ 24 (Romani & Johnston)j PSR 

black body generally peaking m the soft X-ray range, while u-u-j. j- • • c ■ j i ^i 

JO ,i ii ! B1937+21 exhibits sporadic emission of giant pulses m the 

the latter is characterized by a power-law like spectrum , j/cu coi mo en j^i mnf 

J ,,1 , radio band (Sallmcn & Backer 1995: Gognard et al. 1996, 

over a broad energy band. Data available up to now, how- , . , onnn ^ c v. i j. i u j. 

, i i , i i r Kmkhabwala et al. 2000). Such pulses are extremely short 

ever, leave open the debate about the actual presence of , / _ n , n OQ ^, TT n c j . n i 

' . . „ nn . , r ,,rr^ events (r < 0.3 us at 2.38 GHz) confined to small pfiasc 

thermal emission from MSPs. 1 he spectra of the few MSPs . , . .,. . , . . , 

r . , ii. i ii windows trailing the mam and mterpulse. 
lor which thermal models have been claimed, can also be 

X-ray emission from this pulsar was detected by ASCA 

(Takahashi et al. 2001) above 2 keV, with a pulse profile 

Send offprint requests to: G. Cusumano: characterized by a single sharp peak and a pulsed frac- 
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rival times, these authors claimed that the X-ray pulse is 
aligned with the radio interpulse. Later, BeppoSAX de- 
tected pulsed emission from PSR B1937+21 (Nicastro ct 
al. 2002, 2003) and the pulse profile was found to show 
a double peak pattern with a phase separation from PI 
to P2 of 0.52 ± 0.02 and a significance of the second peak 
of ~ 5er. The BeppoSAX data did not allow to study the 
relative alignment between X-ray and radio pulses, be- 
cause the timing did not maintain the necessary accuracy 
to UTC. 

In this letter we present the results of the timing 
and spectral analysis of a RXTE observation of PSR 
B 1937+21. We compare the absolute phases of the X-ray 
and radio pulsed signals and show that the X-ray peaks 
are phase aligned with the radio giant pulses (Section 3). 
The pulsed spectrum is derived in Section 4. 

2. X-ray and Radio Observations and Data 
Reduction 

The RXTE pointings at PSR B1937+21 were performed 
between February 22 and February 27, 2002. The total 
exposure times were about 140,000 s for the PCA units 
0, 2 and 3, and about 20,000 s for the units 1 and 4. 
Standard selection criteria were applied to the observa- 
tion data excluding time intervals corresponding to South 
Atlantic Anomaly passages and when the Earth's limb was 
closer than 10 degrees and the angular distance between 
the source position and the pointing direction of the satel- 
lite was larger than 0?02. We verified that the selection of 
all PCA detector layers, instead of those from the top layer 
only, increased significantly the S /N ratio of the pulsation 
and adopted this choice for our timing analysis. We used 
only data obtained with the PCA (Jahoda et al. 1996) 
accumulated in "Good Xenon" telemetry mode for the 
timing and spectral analysis. Events are time-tagged with 
a 1 (is accuracy with respect to the spacecraft clock and 
with an absolute time accuracy of 5-8 fxs with respect to 
UTC. The UTC arrival times of all selected X-ray events 
were first converted to the Solar System Barycentre us- 
ing the (J2000) pulsar position given in Table 1 and the 
JPL2000 planetary ephemeris (DE200, Standish 1982). 

The radio ephemeris of PSR B1937+21 were obtained 
from high precision timing observations made with the 
100-m Effelsberg radiotelescope in Bonn, Germany, and 
with the Westerbork Synthesis Radio Telescope (WSRT) 
in Westerbork, The Netherlands. From October 1996 tim- 
ing data were collected at 1410 MHz with the Effelsberg 
telescope once per month, with a typical observing time of 
3x7 min. Two circular polarized signals were mixed down 
to an intermediate frequency, detected and coherently 
de-dispersed in a digital filterbank (Backer et al. 1997). 
The data were time stamped with clock information from 
a local hydrogen maser clock and later synchronized to 
UTC(NIST) using the signals from the Global Positioning 
System (GPS). Pulse times-of-arrival (TOAs) were calcu- 
lated from a fit of a synthetic template to the observed 
profile with a template matching procedure (for details 



Table 1. JPL DE200 ephemeris of PSR B1937+21 (data 
from Effelsberg and WSRT observations). 



Parameter 



Value 



Right Ascension (J2000) 

Declination (J2000) 

Proper motion in R.A. (mas yr _1 ) 

Proper motion in Dec (mas yr _1 ) 

Frequency (Hz) 

Frequency derivative (Hz s _1 ) 

Frequency 2nd derivative (Hz s~ ) 

Epoch of the period (MJD) 

DM (cm" 3 pc) 

DM derivative (pc cm~ 3 yr _1 ) 
Epoch validity start/end (MJD) 
Epoch of arrival time' (MJD) 
Frequency of arrival time (MHz) 
Post-fit timing rms (ns) 



19 h 39 m 38!560096(6) 
21° 34' 59'.'13552(13) 
-0.05(2) 
-0.47(4) 

641.928244534462(3) 

-4.331046(13) x 10~ 14 

8.9(1.8) x 10~ 27 

52328.0 

71.02666(16) 

-2.10(17) 

50360 - 52867 

52332.167314839268(11) 

1409.300 

910 



Note: Uncertainties quoted are in the last digit (s) and 
represent 2a estimates (twice the formal tempo errors). The 
Epoch of arrival time refers to the maximum of the main 
peak in the integrated radio profile, 
t at Effelsberg telescope 



see Lange et al. 2001). TOA errors for PSR B1937+21 lie 
in the range 80-200 ns making it one of the most precise 
objects for pulsar timing. 

The WRST observations were performed since July 
1999 at center frequencies of 840 and 1380 MHz using 
the Dutch pulsar machine PuMa (Voute et al. 2002). The 
observing and analysis procedure to produce the WSRT 
TOAs were similar to those at Effelsberg except that lin- 
ear polarisations were recorded. The dual frequency na- 
ture of the WSRT data set allowed us to accurately moni- 
tor dispersion measure (DM) variations which could cause 
significant uncertainties in the absolute timing of pulsars 
(Backer et al. 1993). 

In the timing analysis, both sets of TOAs obtained at 
Effelsberg and WSRT were first independently fit to a pul- 
sar spin-down model with the software package tempo 1 . 
The resultant radio ephemcrides were then used for align- 
ing the RXTE data with the radio profiles, producing fully 
compatible results. Finally, we produced a best-fit timing 
model for PSR B1937+21, as given in Table 1, from the 
combined Effelsberg and WSRT TOAs to align the RXTE 
and radio data (see Fig. 1). 

3. Pulse profile and phase analysis 

RXTE data were searched for pulsed emission by using 
the folding technique in a range centered at a value com- 
puted from the ephemeris reported in Table 1. The plot 
of the x 2 vs the pulsar frequency showed a clear single 
maximum, very prominent above the noise level and the 
pulsar frequency, estimated by fitting the x 2 peak with 
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Fig. 1. Top: The RXTE pulse profile of PSR B1937+21 
in the 2-17 keV energy band. The bin size corresponds 
to - 16 (is. Middle: The 4-10 keV BcppoSAX profile 
(Nicastro ct al. 2002) with the PI aligned with the PI 
phase in the top panel. Bottom: Radio pulse profile at 
1.6 GHz obtained from Effelsberg. Vertical arrows indicate 
the phases of observed radio giant pulses. 



a gaussian profile, was v = 641.92824453 ± 2 x 10~ 8 Hz 
in agreement within the errors with that from the radio 
timing model (Table 1). 

The highest significance of the pulsation is reached in 
the energy interval 2.5-17.0 keV and the resulting X-ray 
pulse profile, obtained by folding the data with the radio 
frequency (Table 1), is shown in Fig. 1 (top panel). It is 
characterized by a prominent narrow first peak (PI), and 
by a less apparent second peak (P2). The significance of 
the latter is 5a above the off-pulse level. P2 is lagging PI 
by 0.526 ± 0.002, determined by fitting both pulses with 
symmetric Lorentzian shapes. The detection of P2 is con- 
firmed by the BcppoSAX data in which a second peak 
was detected at the same phase location (Nicastro et al. 
2002). We show the BeppoSAX result in the middle panel 
of Fig. 1: this profile has been shifted to align PI with the 
phase of PI in the top panel. Pulse widths are wider in 
the BeppoSAX data likely because the events are affected 
by a less accurate time tagging. The PI width (FWHM) 
measured in the RXTE profile is only 29 ± 2 /is and the P2 
width is 51 ± 21 [is. Fig. 1 (bottom panel) shows the ra- 
dio profile from one Effelsberg observation. Vertical arrows 



mark the phases of giant pulses (Kinkhabwala et al. 2000). 
The comparison in absolute phase between the X-ray and 
radio profiles shows that the PI lags the main radio pulse 
by 44 ± 1 ± 5 ± 8 /xs and P2 lags the secondary radio peak 
by 51±3±5±8 /is, where the quoted uncertainties are due 
to statistical error, source position inaccuracies and abso- 
lute timing precision of RXTE, respectively. The X-ray 
peaks appear closely aligned with the phase of the radio 
giant pulses. In addition, the phase separation between 
the X-ray pulses of 0.526 ± 0.002 is more consistent with 
the phase separation between the positions of the giant 
radio pulses (0.5264 ± 0.0006) than with that between the 
radio main and secondary pulses (0.52106 ±0.00003). The 
latter makes a systematic difference in the absolute X-ray 
and radio timing as explanation for the shifts unlikely. 

The occurrence of the same phases for the X-ray pulses 
and the radio giant ones suggests the possibility that high 
energy photons are emitted simultaneous with the lat- 
ters. Therefore, we searched if there is some evidence for a 
bunching of X-ray photons with a rate similar to that of gi- 
ant pulses and equal to ^ 4 pulses per minute (Cognard et 
al. 1996; Kinkhabwala et al. 2000). During the RXTE ex- 
posure we then expected that pulsed events occur in about 
9000 X-ray flares. To investigate this hypothesis we made 
an X-ray light curve selecting only events within the phase 
interval centered in PI with a phase width of A<fi = 0.06 
(90 /is) and studied the frequency distribution of these 
events. Since the dead time of the PCA is about 10 /is 
the maximum content of a bin in the presence of a X-ray 
flares cannot exceed 8-9 counts. We found the following 
statistics: 2 bins with 4 counts, 11 bins with 3 counts, 
574 bins with 2 counts, 294 060 bins with 1 count and 
92 208 884 with count. This distribution is not consis- 
tent with the Poisson statistics, where the expected num- 
ber of bins with a number equal or higher than 2 counts 
is much lower than measured. However, there is no evi- 
dence for the existence of X-ray giant pulses because the 
number of bins with a content different from the Poisson 
distribution was only 116, much lower than the number 
foreseen from the frequency of radio giant pulses. Another 
possibility is that the rate of X-ray giant pulses could be 
lower than that observed in the radio band and that the 
high energy emission could be a mix of steady pulsation 
plus some more rare giant pulse episodes. We constructed 
other light curves selecting events in 10 different phase in- 
tervals far from PI and P2 and with the same phase width 
used in the selection of the PI interval. We found similar 
deviations from the expected Poisson distribution in all 
light curves. In particular, the number of bins deviating 
from a Poisson distribution was found to be between 60 
and 150. Therefore, we conclude that there is no evidence 
that the X-ray emission of PSR B1937+21 is bunched in 
relatively rare events of high intensity. 

4. Spectral analysis 

Pulse-phase histograms were evaluated separately for each 
unit of the PCA for the 256 PHA channels. Response 
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Fig. 2. The pulsed spectrum (PI) with superimposed the 
best fit power law model (top) and their ratio (bottom) 

matrices and sensitive areas were then derived for each 
PCU and summed weighting them by the integrated back- 
ground subtracted counts of the corresponding PCU's 
pulse phase histograms. Pulsed spectra of the main pulse 
were obtained accumulating in the phase interval 0.68- 
0.72 and subtracting the mean off-pulse level evaluated 
in the phase interval 0.27-0.65. Because of the low inten- 
sity of P2, its phase interval was not added to the pulsed 
signal. Spectra were finally combined by summing the in- 
dividual PHA counts and assigned a total exposure time 
equal to the sum of the individual exposures. Furthermore, 
energy channels were rebinncd to have a minimum content 
of 20 counts. The resulting pulsed spectrum in the range 
2-25 keV was fitted with an absorbed power-law model 
with the interstellar hydrogen column density Ah fixed 
to 2.0 x 10 22 cur 2 (Nicastro et al. 2003). The best fit 
pulsed spectrum gave a photon index of 1.14±0.07 with 
a reduced \ 2 of 0.90 (47 d.o.f): it is shown in Fig. 2 (top 
panel) with the residuals in units of data to model ratio 
(bottom panel) . These results are in agreement with those 
obtained from the analysis of the data of a BeppoS AX ob- 
servation (Nicastro et al. 2003). The unabsorbed flux in 
the 2-25 keV band is 6.6 x 10~ 13 erg cm -2 s _1 with the 
corresponding luminosity L x = 8.4 x 10 31 9 (d/3.6 kpc) 2 
erg s _1 , where is the pulsar beam size. 

A fit with a black body model gives kT = 2.6±0.1 keV, 
with a, Xr = 1-5 (47 d.o.f). Although the fit is marginally 
acceptable, there are some systematic deviations in the 
residuals, and the derived temperature T = 2.8 x 10 7 K is 
very high, particularly for an old recycled MSP. 

5. Discussion 

The RXTE observations of the millisecond pulsar PSR 
B 1937+21 provided for the first time a detection of a 
pulsed emission up to an energy of 20 keV. Like in the ra- 
dio, we find a double peaked pulse shape with a dominant 
first peak and a much weaker second peak. Despite being 
only weakly detected, its presence is confirmed by recent 
BcppoSAX data showing a second peak at the exact same 
location. The pulsed spectrum is remarkably hard with a 
photon index of 1.14 indicating that the X-ray emission 



is non-thermal. The major finding of the present analy- 
sis is that the X-ray and radio pulses are not precisely 
aligned, but the formers lag the latters by a phase dif- 
ference of 0.028 ± 0.006 (PI). The X-ray pulses are then 
very well aligned to the giant pulses observed in the radio 
band (Kinkhabwala & Thorsett 2000). Our result is not in 
agreement with the single peak profile at the same phase 
of the radio interpulse, reported by Takahashi et al. (2001) 
from ASCA data. 

Phase coincidences between radio giant and X-ray 
pulses have earlier been reported by Romani & Johnston 
(2001) for the MSP PSR B1821-24. These authors pre- 
dicted that PSR B1937+21 must also show the same prop- 
erty, which has now been confirmed. Among the ordinary 
pulsars, giant pulses are observed only from the Crab pul- 
sar (e.g. Lundgrcn et al. 1995, Cordes et al. 2003) and 
very recently from PSR B0540— 69 (Johnston & Romani 
2003). 

The fact that X-ray and giant radio pulses are observed 
in the same narrow phase intervals suggests that they are 
emitted from the same region of the magnetosphere. It is 
indeed very unlikely that travel time and aberration effects 
will combine themselves to give exactly the same phases. 
It is unknown if this region lies close to the polar caps or 
in some other region of the magnetosphere. According to 
the scenario proposed by Romani & Johnston (2001) these 
pulses are produced in the outer gaps where there is co- 
pious production of secondary e + e~ pairs. In this respect 
it is important to note that PSR B1937+21 is the only 
known pulsar with an estimated magnetic field strength 
at the light cylinder larger than that of the Crab pulsar, 
while PSR B1821-24, has a field at the light cylinder sim- 
ilar to the Crab. Also PSR B0540— 69 ranks near the top 
of the distribution of pulsars with strong B fields near 
the light cylinder (about half the strength of the Crab 
and PSR B1821-24). The X-ray pulsed emission should 
be characterized by a stable intensity, in fact we do not 
find evidence for any clustering in time that could indi- 
cate the presence of X-ray flares. An interesting possi- 
bility is that there are large instabilities affecting mainly 
the space distribution of secondary pairs rather then their 
number. Likely, the current across the gap should remain 
nearly constant. It is then possible that the amplitude of 
these spatial fluctuations may occasionally be very large 
to produce an enhanced coherence in the radio emission 
observed as giant pulses. We expect that these fluctua- 
tions occur over a very short distance scale, comparable 
or smaller than cAfyjv ~ 1.3 x 10 5 cm, equal to a fraction 
of about 4.4xl0~ 4 of the radius of the velocity-of-light 
cylinder, to keep the phase window of giant pulses of the 
order of Air = 2.8 x 10 -3 , corresponding to 1° as found 
by Kinkhabwala & Thorsett (2000). 
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